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O ■ Abstract 

The isovector dipole response in 208 Pb is described in the frame- 
work of a fully self-consistent relativistic random phase approximation. 
The NL3 parameter set for the effective mean-field Lagrangian with 

^J | nonlinear meson self-interaction terms, used in the present calcula- 

tions, reproduces ground state properties as well as the excitation en- 
ergies of giant resonances in nuclei. In addition to the isovector dipole 
k> | resonance in 208 Pb, the present analysis predicts the occurrence of 

^h '_ low- lying El peaks in the energy region between 7 and 11 MeV. In 

particular, a collective state has been identified whose dynamics corre- 
spond to that of a dipole pygmy resonance: the vibration of the excess 
neutrons against the inert core composed of equal number of protons 
and neutrons. 



The multipole response of nuclei with large neutron excess has been the 
subject of many theoretical studies in recent years. Exotic nuclei with ex- 
treme neutron to proton ratio exhibit many interesting and unique structure 
phenomena. In addition to exotic ground state properties, the onset of low- 
energy collective isovector modes has been predicted and experimental evi- 
dence for these modes in light nuclei has been reported [ij . Catara et al. have 
studied the low-lying components in strength distributions of weakly bound 
neutron-rich nuclei 0, and the effect of large neutron excess on the dipole 
response in the region of the giant dipole resonance in O and Ca isotopes J3]. 
They have shown that the neutron excess increases the fragmentation of the 
isovector giant dipole resonance (GDR) and that the radial separation of 
proton and neutron densities leads to non-vanishing isoscalar transition den- 
sities to the GDR states. The fragmentation of the isoscalar and isovector 
monopole strength in neutron rich Ca isotopes has been studied in Ref. [|], 
and in Ref. |J the onset of pygmy dipole resonances in Ca isotopes has been 
analyzed. 

The pygmy dipole resonance, which is also the subject of the present 
study, results from the excess neutrons oscillating out of phase with a core 
composed of equal number of protons and neutrons. A number of theoretical 
models have been applied in studies of the dynamics of pygmy dipole reso- 
nances. These include: the three-fluid hydrodynamical model (the protons, 
the neutrons of the same orbitals as protons, and the excess neutrons) ||, the 
two-fluid (the core fluid and the neutron excess fluid) Steinwedel- Jensen hy- 
drodynamical model 0, density functional theory ||, and the Hartree-Fock 
plus random phase approximation (RPA) with Skyrme forces || ||. More 
recently, large scale shell model calculations have been performed in stud- 
ies of pygmy and dipole states in O isotopes ||, and dipole and spin-dipole 
strength distributions in n Li fT0|| . 

There is also experimental evidence for possible pygmy dipole states in 



208 



Pb. Studies of the low energy spectrum by elastic photon scattering [11 



photoneutron ||12|| , and electron scattering ]T3[ have detected fragmented El 



strength in the energy region between 9 and 11 MeV. The fine structure ex- 
hausts between 3 and 6% of the El sum rule. The relationship between co- 
herent neutron particle-hole (p-h) excitations and the onset of dipole pygmy 
resonances in 208 Pb has been investigated in the Hartree-Fock plus RPA 
model ||]. A concentration of strength has been found around 9 MeV ex- 
hausting 2.4% of the El sum rule. In particular, two pronounced peaks have 
been calculated at 8.7 MeV and 9.5 MeV, which appear as likely candidates 
to be identified as pygmy resonances. The exact location of the calculated 
pygmy states will, of course, depend on the effective nuclear interaction. 
Therefore, it would be important to compare the predictions of various nu- 
clear effective forces with experimental data. 208 Pb is a particularly good 
example, since all nuclear structure models have been tested in the descrip- 



tion of ground and excited state properties of this doubly magic spherical 
nucleus. The pygmy dipole resonance can be directly related to the neu- 
tron excess, and therefore the splitting between the GDR and the pygmy 
resonance represents a measure of the neutron skin. Precise information on 
neutron skin in heavy nuclei is essential for the quantification of the isovector 
channel of effective nuclear forces. 

In the present study the isovector dipole response in 208 Pb is described 
in the framework of a fully self-consistent relativistic random phase ap- 
proximation (RRPA). The same effective Lagrangian generates the Dirac- 
Hartree single-particle spectrum and the residual particle-hole interaction. 
The RRPA response functions with nonlinear meson terms have been derived 



in Refs. |I3J , and applied in studies of isoscalar and isovector giant resonances. 
However, the model configuration space did not include the negative energy 
Dirac states, and therefore model calculations did not reproduce the results 
obtained with the time-dependent relativistic mean-field model |15, 16[| . The 



contribution of pairs formed from occupied positive-energy states and empty 
negative-energy states (in the no — sea approximation), is essential for cur- 
rent conservation and the decoupling of the spurious state [T7| . In addition, 
configurations which involve negative-energy states give an important con- 
tribution to the collectivity of excited states. In a recent study JTJ] we 
have employed a fully self-consistent RRPA, including configuration spaces 
with negative-energy Dirac states, to calculate the isoscalar dipole resonance 
structure in 208 Pb. Two basic isoscalar dipole modes have been identified, 
and the discrepancy between the calculated strength distribution and current 
experimental data has been analyzed. 

In Fig. 1 we display the isovector dipole strength distribution in 208 Pb (left 
panel), and the corresponding transition densities to the two states at 7.29 
MeV and 12.95 MeV (right panel). The calculations have been performed 
within the framework of self-consistent Dirac-Hartree plus relativistic RPA. 
The effective mean-field Lagrangian contains nonlinear meson self-interaction 
terms, and the configuration space includes both particle-hole pairs and pairs 
formed from hole states and negative-energy states. The discrete spectrum 
of RRPA states has been folded with a Lorentzian distribution with a width 
of 0.5 MeV. 



The strength distribution has been calculated with the NL3 |T^] param- 
eter set for the effective mean-field Lagrangian. This force has been exten- 
sively used in the description of a variety of properties of finite nuclei, not 
only those along the valley of /3-stability, but also of exotic nuclei close to the 
particle drip lines. Properties calculated with NL3 indicate that this is prob- 
ably the best effective interaction so far, both for nuclei at and away from 
the line of /3-stability. In particular, in Ref. [TIJ it has been shown that the 
NL3 (K nm = 271.8 MeV) effective interaction provides the best description 
of experimental data on isoscalar giant monopole resonances. The calculated 



energy of the main peak in Fig. 1 E p — 12.95 MeV has to be compared with 
the experimental value of the excitation energy of the isovector giant dipole 
resonance: 13.3 ±0.1 MeV [2(|. In the energy region between 5 and 11 MeV 
two prominent peaks are calculated: at 7.29 MeV and 10.10 MeV. In the 
following we will show that the lower peak can be identified as the pygmy 
dipole resonance. 

The transition densities to the states at 7.29 MeV and at 12.95 MeV are 
displayed in the right panel of Fig. 1. The proton and neutron contributions 
are shown separately; the dotted line denotes the isovector transition density 
and the solid line has been used for the isoscalar transition density. As it 
has been also shown in Ref. 0, although the isoscalar B(E1) to all states 
must identically vanish, the corresponding isoscalar transition densities to 
different states need not to be identically zero. The transition densities for 
the main peak at 12.95 MeV display a radial dependence characteristic for 
the isovector giant dipole resonance: the proton and neutron densities os- 
cillate with opposite phases; the total isovector transition density is much 
larger than the isoscalar component; at large radii they both have a similar 
radial dependence. A very different behavior is observed for the transition 
densities to the state at 7.29 MeV: the proton and neutron densities in the 
interior region are not out of phase; there is almost no contribution from the 
protons in the surface region; the isoscalar transition density dominates over 
the isovector one in the interior; the large neutron component in the surface 
region contributes to the formation of a node in the isoscalar transition den- 
sity. In Ref. || it has been shown that this last effect is also characteristic 
for very neutron-rich systems. The transition densities to the state at 10.10 
MeV display a radial behavior which is intermediate between those shown in 
Fig. 1, but closer to the isovector giant dipole at 12.95 MeV. We have also an- 
alyzed the RRPA amplitudes of the three states: the neutron p-h excitations 
contribute 65%, and the proton 35% to the total intensity of the isovector 
giant dipole at 12.95 MeV, while the neutron contribution is 86% for the 
state at 7.29 MeV. The proton p-h excitations contribute only 14% to the 
total RPA intensity of this state. For the state at 10.10 MeV, on the other 
hand, we find 68% of proton excitations and only 32% is the contribution 
from neutron p-h configurations. However, 31% of the total intensity comes 
from a single proton p-h state: gl '/2~ 1 h9/2. We notice that in the study 
of neutron halos and El resonances in 208 Pb ||, performed in the HF+RPA 
model with the SGII interaction, it was found that for the pygmy states the 
neutron response is a factor 10 larger than the proton response, whereas at 
energies corresponding to the GDR this ratio is about 1.6 or roughly N/Z. 

The phenomenon of low-lying isovector dipole strength was already stud- 
ied almost thirty years ago in the framework of the three-fluid hydrodynam- 
ical model 0. By using a generalization of the Steinwedel- Jensen model [^T] 



to three fluids: the protons, the neutrons in the same orbitals as protons, 



and the excess neutrons, two normal modes of dipole vibrations were identi- 
fied: (i) vibrations of the protons against the two types of neutrons, and (ii) 
the vibration of the excess neutrons against the proton-neutron core. In the 
case of neutron-rich nuclei, the later mode corresponds to pygmy resonances. 
For 208 Pb, in addition to the GDR state at 13.3 MeV, a low-lying pygmy 
state at 4.4 MeV excitation energy was found in the analysis of Ref. ||. The 
dipole strength of this state, however, was negligible (2 orders of magnitude) 
compared to the GDR state. 

In Fig. 2 we plot the transition densities to the two states at 7.29 MeV and 
10.10 MeV. The contributions of the excess neutrons (82 < N < 126) (solid), 
and of the proton- neutron core (Z, N < 82) (dashed) are displayed separately. 
By comparing with the transition densities shown in Fig. 1, we notice that 
there is practically no contribution from the core neutrons (N < 82). The 
reason is, of course, that the p-h configurations which involve core neutrons 
have much higher excitation energies. For the GDR state at 12.95 MeV 
the transition densities of the excess neutrons and the core have the same 
sign in the interior (r < 3/m), and opposite phases in the surface region. 
The overall radial dependence is, however, very similar. The two transition 
densities have the same sign for the state at 7.29 MeV. The core contribution, 
however, vanishes for large r and only oscillations of the excess neutrons are 
observed on the surface of 208 Pb. 

The difference in the collective dynamics of the two modes is also ex- 
emplified in the study of transition currents. In Figs. 3 and 4 we plot the 
velocity fields for the peaks at 7.29 MeV and 12.95 MeV, respectively. The 
velocity distributions are derived from the corresponding transition currents, 
following the procedure described in Ref. [p^] . In both figures the velocity 
field of the proton- neutron core {Z,N < 82) (left panel), is separated from 
the contribution of the excess neutrons (82 < N < 126) (right panel). To the 
largest velocity in Figs. 3 and 4 a vector of unit length is assigned. All the 
other velocity vectors are normalized accordingly. We notice that both the 
core and the excess neutrons contribute to the velocity field of the giant reso- 
nance state (Fig. 4), though the largest velocities correspond to the vibrations 
of the excess neutrons in the surface region. For the state at 7.29 MeV, on 
the other hand, the core velocities are much smaller than those of the excess 
neutrons on the surface. The velocity fields in Fig. 3 corroborate the picture 
of dipole pygmy resonances as oscillations of the excess neutrons against the 
inert core of protons and neutrons in the same shell model orbitals. 

In conclusion, we have used a fully self-consistent relativistic random 
phase approximation to analyze the isovector dipole resonance structure in 
208 Pb. In particular, we have investigated the relationship between the neu- 
tron particle-hole excitations at low energy and the onset of dipole pygmy 
resonances. For the effective interaction we have used the set of nonlinear pa- 
rameters NL3. This interaction has been used in many recent applications of 



the relativistic mean field model, both in stable nuclei and in nuclei far from 
the valley of /3-stability, including drip-line systems. In the particular case 
studied in this Letter, model calculations reproduce the excitation energy of 
the isovector giant dipole resonance in 208 Pb. In addition, low-lying El peaks 
are calculated in the energy region between 7 and 11 MeV . While some of 
these represent the fine structure of the giant resonance, a collective state has 
been identified whose dynamics correspond to that of a dipole pygmy reso- 
nance. By analyzing transition densities and velocity distributions, we have 
related the onset of pygmy resonance at 7.29 MeV to the vibration of the 
excess neutrons against the inert core composed of protons and neutrons in 
the same shell model orbitals. Experimental information on the fragmenta- 
tion of El strength in 208 Pb is, however, only available in the energy window 
8 — 12 MeV. In order to test the predictions of the present analysis, it would 
be important to obtain experimental data also in the energy region below 8 
MeV. 
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Figure Captions 

• Fig.l Isovector dipole strength distribution in 208 Pb (left panel), and 
transition densities for the two peaks at 7.29 MeV and 12.95 MeV (right 
panel). Both isoscalar and isvoector transition densities are displayed, 
as well as the separate proton and neutron contributions. All transition 
densities are multiplied by r 2 . 

• Fig. 2 Isovector dipole transition densities to the 7.29 MeV and 12.95 
MeV RRPA states in 208 Pb. The contributions of the excess neutrons 
(82 < N < 126) (solid), and of the proton- neutron core (Z,N < 82) 
(dashed) are displayed separately. The transition densities are multi- 
plied by r 2 . 



• 



Fig. 3 Velocity distributions for the RRPA state at 7.29 MeV in 208 Pb. 
The velocity field of the proton- neutron core (Z, N < 82) (left panel), is 
separated from the contribution of the excess neutrons (82 < N < 126) 
(right panel). 

Fig. 4 Same as in Fig. 3, but for the RRPA state at 12.95 MeV in 
208 Pb. 
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